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ABSTRACT 

We present Australia Telescope Compact Array (ATCA) and Robert C. Byrd 
Green Bank Telescope (GBT) observations of ammonia (NH 3 ) and the 1.2 cm 
radio continuum toward the ultraluminous infrared galaxy (ULIRG) merger 
Arp 220. We detect the NH 3 (1,1), (2,2), (3,3), (4,4), (5,5), and (6,6) inversion 
lines in absorption against the unresolved, (62 ± 9) mJy continuum source at 
1.2 cm. The peak apparent optical depths of the ammonia lines range from ~ 
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0.05 to 0.18. The absorption lines are well described by single-component Gaus- 
sians with central velocities in between the velocities of the eastern and western 
cores of Arp 220. Therefore, the ammonia likely traces gas that encompasses both 
cores. The absorption depth of the NH 3 (1,1) line is significantly shallower than 
expected based on the depths of the other transitions. The shallow (1,1) profile 
may be caused by contamination from emission by a hypothetical, cold (< 20 K) 
gas layer with an estimated column density of < 2 x 10 14 cm~ 2 . This layer would 
have to be located behind or away from the radio continuum sources to produce 
the contaminating emission. The widths of the ammonia absorption lines are 
~ 120 — 430 km s _1 , in agreement with those of other molecular tracers. We 
cannot confirm the extremely large linewidths of up to ~ 1800 km s -1 previously 
reported for this galaxy. Using all of the ATCA detections except for the shallow 
(1,1) line, we determine a rotational temperature of (124 ± 19) K, correspond- 
ing to a kinetic temperature of T^ n = (186 ± 55) K. Ammonia column densities 
depend on the excitation temperature. For excitation temperatures of 10 K and 
50 K, we estimate N(NH 3 ) = (1.7 ±0.1) x 10 16 cm- 2 and (8.4 ±0.5) x 10 16 cm- 2 , 
respectively. The relation scales linearly for possible higher excitation tempera- 
tures. Our observations are consistent with an ortho-to-para-ammonia ratio of 
unity, implying that the ammonia formation temperature exceeds ~ 30 K. In the 
context of a model with a molecular ring that connects the two nuclei in Arp 220, 
we estimate the H 2 gas density to be ~ fy ' 5 x (1 — 4) x 10 3 , where / v is the 
volume filling factor of the molecular gas. In addition to ammonia, our ATCA 
data show an absorption feature adjacent in frequency to the NH 3 (3,3) line. The 
line does not appear in the GBT spectrum. If we interpret the line to be from the 
OH 2 n 3 / 2 J = 9/2F = 4— 4 transition, it would have a linewidth, systemic veloc- 
ity, and apparent optical depth similar to what we detect in the ammonia lines. 
Comparing the new line to the previously detected 6 GHz OH 2 n 3 / 2 J = 5/2 
F = 2 — 2 transition, we determine a rotational OH temperature of ~ 245 K, 
about 2 times the rotational temperature of ammonia. If this association with 
OH is correct, it marks the first detection of the highly excited (~ 511 K above 
ground state) 2 n 3 / 2 J = 9/2 F = 4 — 4 OH line in an extragalactic object. 

Subject headings: galaxies: individual (Arp 220) — galaxies: starburst — galax- 
ies: ISM — galaxies: nuclei — ISM: molecules — radio lines: galaxies 



Arp 220 
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Introduction 



Ultraluminous infrared galaxies (ULIRGs), with far-infrared luminosities exceeding 
-^fir <i 1O 12 L , exhibit so me of the highest star fo rmation (SF) rates encountered in the 
Universe (for a review, see ISanders fc Mirabellll996l ). They are mergers containing at least 
two massive, gas-rich galaxies whose dust is heated during the mergin g process. Cloud- 
cloud collisions and the subsequent inflow of gas feed nuclear activity ( jNagar et al.l 120031 ) 
and trigger an amazingly high rate of star f ormation in the c entral few 100 pc, typically 



reaching tens to hundreds of M yr 1 (e.g., iFlores et al.l 120041 ) . ULIRGs are particularly 



important galaxies as they are the most readily detectable p opulation for which the SF per 



co-moving volume can be traced out to high redshifts (e.g., iMadau et al.l 11998c I Cole et al. 



200 ll ; iDarling &: Giovanellil 120021 ) . Furthermore, the number of luminous IR galaxies per 



p~ - „ , _ 

Bridge et al. 


2007; 


Hopkins et al. 


2010) 



duction of stars early in the Universe (e.g. 
increasing evidence that U LIRGs may enter an optically bright quasar phase subsequently in 



the merging process (e.g. iHutchings fc Neflj Il988l ; ICanalizo Stockton! 1200 ll ; iDasyra et al. 
2007c IVeilleux et al-lEooa ). 



The prototypical ULIRG is Arp220 (IC4553; L F i R ~ 1.5 x 1O 12 L ), at a distance 
of only D ~ 75 Mpc. Its proximity allows us to study physical processes in more detail 



than we can for its more dista nt cousins. Arp220 has two starburst nuclei (e.g., iNorris 



1988c iDownes &; Solomonlll998l ) separated by / . / 95 (350 pc) at a position angle of 92°. Many 



-™ „~ . „^ „ — „. „, „ „ , ^, ^, ^, 

HCN, HNC, HCO+, H 9 C0, H 9 0,H 9 S, HC,N, HNCO, and thermal OHfe.fr.. Solomon et al. 


1990l. 1992: Hiittemeister et al. 1995a: IScoville. Yun. & Brvant 


1997; ISolomon et al. 


1997; 


Gerin & Phillips! 


1998[ 


PaDadopoulos & Seaouistl 19981: Yao et al. 


2003; 


Arava et al. 


2004; 


Downes & Eckart 


2007 


ISalter et al. 20081: Aalto et al. 20091: 


Greve et al. 


2009: 


Vlartm et al. 


2011) as well as OH megamaser emission (e.g., Baan et al. 


198: 


>; 


Rovilos et al 


. 12003) 


. The 



molecular gas traced by (sub-)arcsecond CO imaging reveals small, dense disks around each 
of the two nuclei. These small disks are embedded in a ring or disk that encompasses both 
cores, and that ring is itself embe dded in an even larger molecular disk, ~ 2.5 kpc in size 



and with a density of ~ 10 2 cm 3 (|Scovi 



Sakamoto et al. 1999; Mundell et al. 



le. Yun. fc Bryantlll997l ; IDownes fc Solomon! 11998 



20011 ) . The western nucleus has a hig her molecular mass 



than its eastern counterpart, and may contain a supermassive black hole (IDownes fc Eckart 



20071 ). Both nuclei show signatures of outflows (jSakamoto et al.ll2009l ). Direct measurements 
of the physical properties of the gas, however, are difficult. For linear molecules, radiative 
transfer models, like the Large Velocity Gradient (LVG) analysis, have solutions that are 
degenerate between high temperature/low density and low temperature/high density. Anal- 
ysis of the excitation conditions show that the two cores have temperatures in the range 
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45 — 120 K and densi ties of order 10 6 cm 3 , with the we s tern nucleus being som ewhat 



warmer and denser (e.g. lAalto et al.ll2009t iGreve et al.ll2009t iMatsushita et al.ll2009l ). Far- 



ther out, while the gas is known to be warm (T > 40 K), no accurate temperature or density 
measurements exist. 



Solomon et al.l ( 119921 ) and iGao fc Solomonl (l2004ai lbl) observed CO and HCN in a large 
sample of ULIRGs. They conclude that HCN is a much better tracer of SF than CO, and 
in fact they find that the HCN line luminosity scales linearly with the SF rate in galaxies. 
HCN exhibits a large electric dipole moment and the molecule is only excited when the gas 
is sufficiently dense (n(Hz) > 10 4 cm -3 ). Dense molecular g as is therefore the best tr acer 
for the material that is readily converted into stars (but see iGracia-Carpio et al.l 120061 for 
caveats). 

Similar to HCN, ammonia (NH 3 ) also traces dense gas, though it is excited at densities 
down to ~ 10 3 ~ 4 c m" 3 . Its abundance is closely linked to that of HCN by the cyanide chem- 
istry network (e.g., ISchilke et al.lll992l ). In contrast to linear molecules like CO and HCN, the 
ammonia molecule is a prolate symmetric top with tetrahedral structure. Most of its excita- 
tion states, characterized by the quantum numbers (J,K), show 'inversion doublets' caused 
by the nitrogen atom tunneling through the plane defined by the three hydrogen atoms. The 
lowest inversion doublet of each .ff -ladder (with the quantum numbers J—K) is metastable. 
The relative populations of the metastable doublets are mainly determined by collisions, 
follow a Boltzmann distribution, and thus provide one of the most suitable thermometers to 
estimate kinetic temperatures of extragalactic mol e cular gas (e.g., IWalmslev &: Ungerechts 



19831 : lUngerechts. Winnewisser. fc Walmslevl Il986l : banbv et al.l Il988l : iHenkel et all boooh . 



Since the energy separation between the individual inversion doublets is only a weak func- 
tion of J and K, a vast range of excitation properties can be covered by observations of 
lines closely spaced in frequency. The most relevant NH 3 lines fall in the 1.2 cm band 
[y ~ 25 GHz). Relative calibration of these lines tends to be very good since they can be 
observed nearly simultaneously, under similar atmospheric conditions, and with the same 
telescope, receiver, and backend. 



NH 3 has been detected in Arp220 by lTakano et al.l (120051 ) with the 45 m Nobeyama sin- 
gle dish telescope. Their main result was that the ammonia has linewidths of ~ 1800 km s _1 , 
which they suggest could indicate that the gas is located in an accretion disk around an ac- 
tive nucleus. Alternatively, the wide lines may indicate strong molecular inflows or outflows. 
Single dish observations, however, can be affected by instabilities that lead to poor spectral 
baselines and uncertain line profiles. This is particularly the case for broad lines. Here, we 
present ATCA interferometric data with higher signal to noise and better spectral baselines. 
We also present a sensitive GBT spectrum. In Sect. [2] we describe our observational setups. 
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W e present our results in Sect.0 along with a critical comparison of our data with those 
of iTakano et al.l ( 120051 ). In Sect. 0] we derive and interpret the apparent optical depths, gas 
temperatures, gas densities and kinematics of Arp220, as well as its SF properties. Sect. [5] 
summarizes our paper. 



2. Observations and Data Reduction 

2.1. ATCA Observations 

We observed the spectral lines of ammonia and the 1.2 cm (K-band) continuum toward 
Arp 220 with the ATCA0. The data were taken in several stages: (1) A broadband spectrum 
was obtained covering the 23.2-23.6 GHz frequency range. At the time of the observations 
(2005 July), the ATCA could not cover this full frequency range in a single observation. 
Therefore, we configured the observations with 128 MHz bands of 64 channels each, and 
stepped through the desired frequency range in 60 MHz increments (indicated by the bars at 
the bottom of Fig.[TJ), with two bands observed at a time. The observations were performed 
over 4 days with ~ 4.5 h on-source per day. The spectral resolution is equivalent to 2.2 
channels (~ 57 km s _1 ). Calibration was performed individually for each frequency band and 
observing day. At the time, the ATCA was in its most compact H 75 antenna configuration, 
where baselines reach from 31 m to 89 m. (2) The (6,6) line was observed in 2005 March with 
the H214 antenna configuration, which has 82m-247m baselines. The integration times on 
source were ~ 2.5 h (2005 March 17), and ~ 4.5 h (2005 March 20). (3) On 2005 July 30 
we observed the (4,4) line on source for 4.5 h in the H75 configuration. (4) A day later we 
measured the (5,5) line with the same integration time in the same array configuration. The 
correlator settings for stages (2)-(4) were the same as in stage (1), except that, for the (4,4) 
and (5,5) transitions, only two 128 MHz wide frequency settings were used in parallel and 
separated by 72 MHz. The (6,6) line was observed in a single frequency window. A summary 
of the observational setups is provided in TableHJ 

The data were processed using the MIRIAD software package. We observed PKS 1253- 
055 for bandpass calibration, PKS 1934-638 for flux calibration, and PKS 1538+149 (in 
20min intervals) for phase/gain calibration. For each spectrum, we subtracted the contin- 
uum as determined by a linear fit to the line-free channels. In the process we also created 
continuum data sets. All data were Fourier transformed using 'natural' weighting. To obtain 



1 The Australia Telescope Compact Array is part of the Australia Telescope which is funded by the 
Commonwealth of Australia for operation as a National Facility managed by CSIRO. 
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the spectra, the data were not deconvolved due to the lack of spatial structure on the scales of 
the dirty beam (source size -C 10" versus ~ 10" — 40" beam size, the latter depending on the 
array configuration, see Tabled]). Deconvolution with Clean is only needed to reconstruct 
source morphologies. This, however, comes for the prize of introducing non-linearities. For 
our unresolved source, we decided not to Clean the data. Given the relatively strong con- 
tinuum flux, this implies that the noise within the planes of the dirty images is dominated by 
sidelobes and we measure 5 — lOmJybeam -1 depending on the observations. For an estimate 
of the more representative statistical noise and thus for the apparent optical depth uncer- 
tainty, however, we did apply Clean deconvolution to remove the sidelobes. This method 
yields rms values of ~ 0.4 — 1.2 mJy beam -1 per channel for the different observations. We 
estimate the absolute calibration uncertainties to be of order 15%. 



2.2. GBT Observations 

We also observed the 1.2 cm ammonia spectrum toward Arp 220 on 2006 January 21 and 
2006 February 02 with the Green Bank Telescope (GBT) of the National Radio Astronomy 
Observatory^ (see Table[T]). We observed both circular polarizations simultaneously using a 
total-power nodding mode in which the target position was placed alternately every 2.5 min 
in one of the two beams of the 22 — 26 GHz K-band receiver. The beams are spaced by 3' 
in the cross-elevation direction. Total time on source was 40 minutes. We configured the 
spectrometer with an 800 MHz spectral window centered on the NHs(3,3) line for Arp 220. 
The data were calibrated and averaged in GBTIDL. During calibration, we used an estimate 
of the atmospheric opacity obtained from weather conditions. We boxcar averaged the 
spectrum to achieve 1.56 MHz channel spacing, and subtracted a fifth-order polynomial fit 
to the line-free channels over the full 800 MHz, a conservative approach for broadband GBT 
data. Finally we combined the spectra from the two days and the rms noise of the final data 
is 0.8 mJy per 1.56 MHz channel. 



2 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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3. Results 

3.1. 1.2 cm Continuum Emission 

With the ATCA, we detect an unresolved continuum source toward Arp 220 at RA ( J2000) 
= 15 h 34 m 57.2 s , DEC(J2000) = 23° 30' 11". The source has a flux density at 23.4 GHz of 
208 mK or 62 mJy (error of order 15%) measured within the 31.1" x 21.1" beam. From March 
to August 2005, the measured continuum flux varied by less than 5%. 



3.2. Ammonia Inversion Lines 

In Figs.[T]and[2]we show ATCA spectra of Arp 220 in units of apparent optical depth (see 
Sect. l4.lj ). and Fig.[3] shows a GBT spectrum overlaid on the ATCA data. All six observed 
ammonia inversion lines, NH 3 (J,K) = (1,1), (2,2), (3,3), (4,4), (5,5), and (6,6), are detected 
in absorption against the radio continuum of Arp 220. We do not detect clear emission 
features, but the broadband ATCA spectrum displays slightly positive values adjacent to 
the (1,1) absorption line (see Sects. l4~TJ and 14.21 for a discussion). We fitted single- component 
Gaussian profiles to all observed absorption lines. The resulting apparent optical depths 
t (see Sec. 14.11 for their derivation), the optical, barycentric peak velocities v^ haTy , the 
linewidths Avi/ 2 , and the apparent optical depths integrated over the line profiles, f rdv, 
are given in Tabled and the fits are shown in Fig[21 The fitted line centroids are all within 
~ 60km s _1 of the systemic velocity of 5434km s -1 (z = 0.018126). This is remarkable 
given the ~ 57km s _1 velocity resolution of the ATCA observations (see Sect. [2]). The 
linewidths range between ~ 230 km s _1 and ~ 330 km s" 1 for most lines. The (1,1) and 
(3,3) lines, however, have widths of 121 ± 38 km s^ 1 and 435 ± 19 km s _1 , more extreme 
than the other transitions. Note, though, that the (3,3) line slightly overlaps with another 
absorption feature, described in Sect. l3.3[ Our fitted lin ewidths are in agreement with those 



fo und for other molecu lar tracers (e.g. lSalter et al.ll2008l ) but contradict the results presented 



m 



Takano et al.l (120051 ) for the NH 3 (1,1) and (3, 3) lines. We cannot c onfirm the extremely 



wide lines of up to ~ 1800 km s 1 observed by iTakano et al.l (120051 ) with the Nobeyama 
45 m telescope. None of the lines measured with the ATCA or GBT have widths exceeding 
~ 450km s" 1 (see Tableland Figs-HHS]). 

There are some notable differences in widths and strengths of the lines measured with 
the GBT and the ATCA. The reason may be the weakness of the lines, the complexity 
of the spectra, and potential confusion with the continuum emission from Arp220. The 
NHs(l,l) and (2,2) lines are separated by less than 30 MHz so, for the GBT data, it is 
difficult to precisely define a spectral baseline near these features. The quality of spectral 
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baselines from the interferometric ATCA data, on the other hand, is very good. We would 
like to add, however, that regridding the 128 MHz blocks of our ATCA observations into 
a single spectrum conveys some uncertainty. F or the relatively isolated NH 3 (4,4) line, our 
results agree well with those of iTakano et al.l ( 120051 ). The centroids of the Gaussian line 
fits between the GBT and ATCA also differ by a few tens of km s _1 (Table|2]). This may 
be caused by the limited signal-to-noise ratios of the spectra, despite the relatively small 
errors of the Gaussian fits. Gaussian fitting may have small systematic errors caused by 
the uncertain baseline subtraction and blended lines. The linewidths and peak apparent 
optical depths of the NH3(1,1) line agree between the ATCA and GBT observations. The 
other lines, however, disagree to some extent in linewidth and peak apparent optical depth. 
The (3,3) line is about twice as broad in the ATCA data as in the GBT data and the peak 
apparent optical depths in the GBT data for the NH 3 (2,2) and (3,3) lines are about 1/3 less 
than those in the ATCA spectrum. 



3.3. Potential Detection of OH 2 n 3/2 J = 9/2 

In addition to the (1,1) through (6,6) ammonia absorption lines, the broadband ATCA 
spectrum exhibits an additional absorption feature at ~ 23.39 GHz (Fig. [I]). As with the am- 
monia lines, this additional spectral line appears in both spectra that were taken at that fre- 
quency (the horizontal bars in Fig.[TJshow the tunings of the individual observations). Note, 
however, that the two underlying individual spectra were observed simultaneously which 
implies that they are subject to a common systematic uncertainty and are not entirely inde- 
pendent measurements. At the redshift of Arp 220, the most obvious line candidates would 
be the F=4-4 and _F=5-5 doublet of the OH 2 n 3 / 2 J = 9/2 transition at rest frequencies of 
23817.6153 MHz and 23826.6211MHz, respect ivehfl This doublet is at an energy ~ 511 K 
above ground state and indicates highly excited gas. Adopting the rest frequency of the 
F=4-4 transition and fitting a Gaussian to this line, we get a line center of —13 ± 10 km s -1 
(relative to the systemic velocity of 5434 km s _1 ) and a linewidth of 292 ± 29 km s _1 , very 
similar to the ammonia lines. These fitting results are also shown in Table|5J Fitting the 
F=5-5 line alone, we find the centroid of the Gaussian would be offset by ~ +100 km s _1 
from systemic, likely too large for a feature that belongs to Arp 220. Since the doublet lines 
are spaced by only ~ 113 km s _1 , however, they may be blended to produce the observed 
line width. A fit from blended lines would have a mean velocity offset of ~ +44 km s _1 - well 
within the range of the line velocities observed. In either case, if the measured absorption 
feature indeed corresponds to one or both OH 2 n 3/2 J = 9/2 lines, it would be the first 



3 taken from the LOVAS catalog: http://physics.nist.gov/cgi-bin/micro/table5/start.pl 
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detection of this transition in an extragalactic object. To date, this doubl et has only been 
observed toward Galactic ultracompact H II regio ns such as W3(OH) (e.g. 
19781 : Iwilson et allll99ll : baudrv fc Mentenlll995h . 



Winnberg et al. 



Lower energy OH doublets have been observed previously in Arp220. These include 

2T 



111/9 J - i- 
(Henkel et al. 


1986. 


1987; 


Salter et al. 


- °l 
2008) 



-20 km s 1 to +5 km s 1 , linewidths of ~ 250—480 km s 1 , and apparent optical depths 
0.05 — 0.11. The parameters of these lines are cons istent with those of our possible 



range 
of r ^ 

ATCA OH line detection. Using our possible line and the ISalter et al.l 120081 2 n 3 / 2 J = 5/2 
F=2-2 detection, we estimate a rotational temperature of ~ 245 K, assuming Local Thermal 
Equilibrium. This estimate is about a factor of two higher than the rotational temperature we 
determine for ammonia (Sect. l4"T2"j) . The line could be verified with wideband interferometric 
observations, for example with the new CABB backend of the ATCA or WIDAR on the 
EVLA. Indirect evidence could come fro m HERSCHEL data, v ia observations of rotational 



lines between different OH doublets (e.g. IWampfler et al.ll2010i ). 



4. Discussion 



4.1. Apparent Optical Depths 



For a molecular cloud that partly covers a continuum source, in the Rayleigh- Jeans 
limit, we can express the measured brightness temperature of a spectral line (Tl) as: 



T L = (fd T cx - f f c T c - f cl T C mb) x (1 - e" 



T 



(IRohlfs fc Wilson! 120041 ) where T ex is the excitation temperature of the line, Tc is the bright- 
ness temperature of a background continuum source, Tcmb is the brightness temperature of 
the cosmic microwave background, and r is the optical depth of the line. The temperatures 
are scaled by the beam filling factors of the cloud (/ c i) and the continuum source (fc)- The 
source covering factor f describes the fraction of the continuum covered by the cloud along 
the line of sight. We use this equation to determine the apparent optical depth r. The source 
covering factor / can be assumed to be close to unity since the individ ual radio continuum 



emitting regions are relatively small, at most a few tens of pc in size (jRovilos et al.l 12003 



Lonsdale et al.ll2006l ; iParra et al.ll2007l ). and supposedly contain mostly supernova remnants 
embedded in the molecular gas (see Sec t. l4.4p . In addition, the mean gas column density in 
the model of iDownes fc Solomon! (119981 ) is Nh 2 ~ 2.5 x 10 24 cm -2 and the mean gas volume 
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density is ~ 500 cm -3 (using the radius of the outer molecular disk of 1.2 kpc, which cor- 
responds for a cylindrical, edge-on geometry to a mean width and line of sight of 1.6 kpc, 
a height of the cylinder of 80 pc, and a gas mass of 5 x 10 9 M Q ; but see Sect. l4.3p . Those 
large densities indicate that molecular gas is abundant and, closer to the starburst cores, 
ammonia is expected to be present over large regions which justifies our assumption of a 
source covering factor near unity. Under these circumstances we can also assume that both 
the cloud and the continuum filling factors are the same for absorption lines, i.e., f c \ = fc- 
Ammonia not along the line of sight toward the radio continuum would appear in emission, 
which is not directly observed (but note the weak, tentative emission features adjacent to 
the NH 3 (1,1) line [Fig.[T] - the peculiarity of the (1,1) transition is discussed in Sect. [472]). 
Under these conditions, Eq.[T] simplifies to 



HTe: 



Tn-T, 



CMBJ 



1 



(2) 



VLBI continuum measurements at 18 cm show that the radio continuum of Arp 220 has 
an extent of < 0.6D" (jSmith et al.lll998l ; iRovilos et al.l 120031 ). which corresponds to a beam 
filling factor of fc < 10~ 3 when compared to the beam of our ATCA broadband observations 
(~ 31'.'1 x 21'.'1). For this value, the measured Arp 220 1.2cm continuum brightness temper- 
ature of T' c = fcTc = 208 mK then implies a true brightness temperature Tq > 208 K. The 
best spatial resolution was obtained for the NH 3 (6,6) line (Sect. l2~Tj) and even for the ~ 8 
times (in area) smaller beam of that observation we do not detect any extent of the con- 
tinuum source. For this smaller beam of the NH 3 (6,6) observation, the same approximation 
then boosts the true continuum temperature limit by the same factor of 8. On the other 
hand, VLBI does filter out flux on larger scales which somewhat counters the low beam 
filling argument. Nevertheless, the true continuum brightness temperature is very likely 
much larger than the CMB temperature of T C mb = 2.73 K. We can therefore neglect Tcmb 
in all further calculations. We can a lso disregard T ex . LVG models of the ammonia molecule 
(e.g.. IWalmsley fc Ungerechtslll983l ) show that, under normal circumstances, the excitation 
of inversion lines hardly exceeds T ex ~ 50 K (see also Sect l4.3|) . The continuum brightness 
temperature therefore dominates the radiative transfer. With those assumptions, Eq.[2]can 
be expressed as 



In 1 - 



\Tl\ 

rpl 

with the approximated r denoting the 'apparent optical depth'. 



(3) 



In Figs.lU [21 and |3] the absorption profiles of the ATCA and GBT observations are 
displayed in units of apparent optical depth (using for both the ATCA continuum flux 
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described in Sec. l3.ip . The Gaussian fits to the spectra (Fig. [2]) result in peak apparent 
optical depths of ~ 0.05 — 0.18 (Table[2]) w hich indicate o ptical ly thin lines. The peak 
apparent optical depths are similar to what iTakano et al.l (120051 ) finds for the NH 3 (1,1), 
(3,3), and (4,4) lines. Our fit to the (2,2) line, however, is about 3 times larger than theirs. 
We should also note that th e integrated apparen t optical depths are different due to the 
larger linewidths reported by lTakano et al.l ( 120051 ). 



4.2. Rotational and Kinetic Temperature 

Using the apparent optical depths, it is possible to derive the ratios of the ammonia 
column densities N(J, K) and excitation temperatures T ex for the two states of a given 
inversion doublet with the equation 

m«l = !.61 x 10» x ££ti> x T x A„ 1/2 (4) 



( IHuttemeister et al.lll995bl ; v in GHz; Avi/2- full width half maximum [FWHM] linewidth 
in kms -1 ). Normalizing these with respect to their statistical weights and assuming that the 
excitation temperatures T ex are the same for the six observed transitions, the ratio of column 
densities and excitation temperatures N/T cx as a function of energy may be described by a 
Boltzmann law with a specific rotational temperature T rot jj/ (or simply Tjj/) via 

W, J')/TL _ 9o P (J') 2J' + 1 exp (z^E\ (5) 



N(J,J)/T ex g op (J) 2J + 1 "V ^JJ' . 
for a pair of inversion lines NH3(J',J') and NH3(J,J) [g op (J) is the statistical weight for a 
given ammonia species, g op = 1 for para-ammonia, g op = 2 for ortho-ammonia]. 

The assumption of equal T ex « T' ex values is justified for t he para-ammonia (1,1), 



\ —1— ; i 

(1983 


) and 


Danbv et al. 


(1988) 



excitation temperatures may vary within a factor of ~ 2 as compared to the lower energy 
para-NH-3 transitions. 

In Fig. HJ we show a Boltzmann plot of the ammonia inversion lines toward Arp 220. The 
GBT and the ATCA data align reasonably well for the NH 3 (1,1) and (2,2) transitions. The 
NH_3(3,3) datum, however, is significantly lower for the GBT observations. The weighted 
population of the NH 3 (2,2) state is larger than that of the (1,1) state which is not expected 
but has been confirmed by our two independent measurements. The small NH 3 (1,1) column 
density is even more surprising as both transitions, NH 3 (1,1) and (2,2), belong to the same 
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para-ammonia species. Both transitions have likely similar excitation temperatures and are 
optically thin (Tabled see also Sect. l4.lj) . This excludes the possibility that saturation effects 
are responsible for different strengths of the (1,1) and (2,2) transitions. 

An explanation for the weak NH3(1,1) feature could be that an additional, cool ammonia 
component in Arp220 contaminates the absorption spectrum by ammonia (1,1) emission. 
This cloud needs to be cool to avoid significant population of the higher energy levels; 
NH 3 (2,2) through (6,6) appear to have absorption profiles that are well described by a single 
rotational temperature and are likely little or not affected by any attenuation from cold 
gas (see the discussion in the last paragraph of this Section). To be seen in emission, the 
cold component would have to be in the beam but away from the line of sight toward the 
continuum. The faint emission channels adjacent to the (1,1) absorption profile of the ATCA 
data may provide an indication for such a cold component. 

In spite of this complication and ignoring the unusual profile of the (1,1) line with 
respect to all other lines, we can estimate the rotational temperature of the gas. Using 
the Boltzmann diagram (Fig.HJ), this is performed by a weighted fit through all (2,2), (3,3), 
(4,4), (5,5), and (6,6) data points (T rot = T 2 34 56 , with an ortho-to-para ammonia abundance 
ratio of 1). Using the ATCA data, we derive T23456 = 124 ± 19 K and the best fit is shown 
in the Boltzmann plot in Fig.Hl Note that the fit describes all data well except for the 
(1,1) ammonia transition. The good fit to the (2,2) though (6,6) data provides confidence 
that our interpretation and Gaussian fits to the related spectra are reasonable in terms of 
peak apparent optical depths and linewidths. However, it also points toward a less reliable 
NH3(3,3) GBT flux measurement. The three lowest transitions of the GBT spectrum do 
not line up at all in the Boltzmann plot. In addition, the GBT NH_3(3,3) datum strongly 
deviates from the best fit of the ATCA NH 3 (2,2) through (6,6) data points. We therefore use 
the superior number and quality of the ATCA spectra throughout the rest of our analysis. 

Since the lowest energetic ammonia (0,0) transition belongs to the ortho-ammonia 
species, cold ammonia formation temperatures would be reflected in an enhanced production 
of ortho-NH 3 over para-NH 3 . In that case one would observe the (3,3) and (6,6) transition 
data to be significantly elevated from the rotational temperature fit through all points. A T ex 
difference between the ortho and para ammonia inversion lines, however, could compensate 
or strengthen a deviation from a single-line fit in the Boltzmann diagram. Such a difference 
can be well fitted by a linear function through NH_3(2,2) and higher, which suggests that the 
ortho-to-para-NH 3 ratio does not deviate substan tially from unity. In turn, this i mplie s that 



the ammonia was formed under warm conditions. [Takano. Nakai. fc Kawaguchil (120021 ) show 



that ortho-to-para-NH-3 ratios near unity require formation temperatures exceeding ~ 30 K. 



For low temperatures, rotational temperatures correspond well to kinetic temperatures. 
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At higher temperatures, however, populations of more and more non-metastable levels as 
well as ra diative decay result in a deviation from the approximation I Walmsley &: Ungerechts 
(cf. 119831 ). In Fig.Owe show the results on the Tk in -T rot relation for different pairs of inver- 
sion transitions . This calculation is based o n an LVG radiativ e transfer model (described in 
Ott et alJliooil using the collision rates of lDanbv et~allll988[ ). For the NH 3 (1,1) and (2,2) 
transitions, T^ n pa T i2 is valid below ~ 40 K. Using higher transitions, the approximate equal- 
ity is valid into progressively higher temperature regimes. For that reason, it is preferable to 
observe high NH 3 transitions if they are bright enough. Nevertheless, even if the T-^ in pa T rot 
relation is broken, one can still apply th e models displa yed in Fig.[5]as they only weakly de- 
pend on volume density (for details, see lOtt et al.ll2005l ). Here we present logarithmic fits of 
Tki n to the rotational temperatures Tyi and T 3 6, that are based on the lowest para- ammonia 
(1,1) and (2,2), and the lowest ortho-ammonia (3,3) and (6,6) transitions, respectively: 



T kin = 6.06 x exp(0.061 T 12 ). 



(6) 



and 



T JT 36 forT 36 <50K 
k " [28.9 x exp(0.015 T 36 ) for T 36 > 50 K 



(7) 



Both fits to the LVG solutions are shown in Fig. [5] as dashed lines to their corresponding 
curves. 

For Arp 220, we determine the rotational temperature across five transitions, rather than 
only a pair. Eq.[6] cannot be applied as the (1,1) level seems to be an outlier in Arp 220 and 
was excluded from our T rot 23 456 fit. We are using Eq. [7] instead. This rotational temperature 
is calculated over a much larger energy difference of the (3,3) and (6,6) levels, rather than 
the smaller difference between (1,1) and (2,2). In addition, a Boltzmann fit through (3,3) 
and (6,6) provides the same 124 K rotational temperature as the fit through (2,2)-(6,6). For 
that reason, we expect that the conversion of Eq. [7J is a good approximation for the kinetic 
temperature of the ammonia-traced gas in Arp 220 and we derive T^ in = (186 ± 55) K. 

Such a kinetic temperature is similar to that fou nd in the more n earby 'normal' starburst 
galaxy NGC253, which exhibi ts T v - m = 130- 210 K Jott et al.ll2005h. Other starburst galax- 



ies, such as Maffei 2, and M 82 (|Hejikel_etLaJ l20o3; Iweifi et al 



2001 



Takano. Nakai. fc Kawaguchi 



20021 ; iMauersberger et al.ll2003l ) show lower rotational temperatures than NGC253, which 
must also result in lower kinetic temperatures. 



For comparison, iDownes fc Eckartl (120071 ) report an intrinsic hot dust temperature com- 
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ponent with Td us t = 170 K toward the western nucleus of Arp220. This is comparable to 
what we obtain for NH 3 , but might not emerge from the same physical feature (see Sect. 14. 4j ). 
Dust temperatures are somewhat lower in the outer molecular ring and over the entire body 
of Arp220. ISO o bservations show two dust temperature components with 47 and 120 K 
(IKlaas et all 119971 ) . The gas thus has a kinetic temper ature component as high as that 



of th e hottest dust component in A rp 220. Recent S MA (ISakamoto et al. 



2009), ASTE fllmanishi et al.ll2010h . and JCMT data (IGreve et al. 



2009 



2008 



Aalto et al. 



Papadopoulos et al. 



2010a() of highly excited, dense gas show typical temperatures of ~ 45 — 120 K in the dens- 
est, ~ 10 6 cm -3 components, likely the gas closest to the cores. The temperature of the gas 
phase we determine is somewhat higher and may arise within the surrounding ring or disk 
(see Sects. I4T31 and l4"~4"t see also Fig.[6]). No direct temperature measurement has so far been 
reported for this physical feature. 

Given the rotational temperatures, we can now estimate properties of any ammonia 
that contributes emission to explain the low absorbed flux of the (1,1) line relative to the 
(2,2) line. Following the Boltzmann plot (Fig.HJ), and assuming that the derived T23456 
rotational temperature can be extrapolated to a hypothetical (1,1) column density, one 
would expect a ~ 4 times larger column density of the (1,1) line than actually measured in 
absorption. According to Eq.|U this would result in an integrated apparent optical depth 
~ 4 times larger. This missing integrated apparent optical depth must be compensated for 
by emission of order ~ 20 mJy. The gas responsible for such emission would have to be cold 
to avoid contaminating the other inversion profiles. For an upper limit, let us assume that 
the NH 3 (2,2) line is uncontaminated at a level two times below the one actually measured 
and that, from that level, the additional emission component raises the (negative) flux to 
the measured value. The emission line would then require a strength of ~ 6 mJy. Using 
this value as an upper limit to the NH 3 (2,2) emission component, we can derive an upper 
limit to the temperatu re from th e (1,1) and (2,2) lines. Using the equations for ammonia 
emission as provided in lOtt et all ( 120051 ). we estimate an upper rotational temperature limit 
of T12 < 20 K. With Eq.[6l this translates to an upper limit of the kinetic temperature of 
the same order of T^ n < 20 K. This scenario can be proven by observations at high spatial 
resolution (< 0.5") that allow direct imaging of emission components offset to the lines of 
sight toward the bright continuum sources. The EVLA may provide such observations in 
its extended configurations. Alternatively, mm/submm interferometers such as ALMA and 
PdBI may be able to confirm the existence of this cold component. The ammonia inversion 
lines, however, do not fall in their wavebands and the predicted cold gas component would 
need to be revealed via low energy, high frequency transitions of other suitable molecular 
tracers. 
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4.3. Gas Densities 



The excitation temperatures of ammonia in ULIRG s is still unclear. An ind i cation of a 
reasonable range might be given by the calculations of IWalmsley fc Ungerechts Jl983h 



sec 

their fig. 3, and using their own collision rates). For the ratio of an NH 3 concentration to 
velocity gradient of 10 _4 cm -3 (km s _1 pc -1 ) -1 , they show that the excitation temperatures 
of the ammonia inversion lines span T cx ~ 10 — 50 K, except for a poorly confined range 
of parameters where the ammonia (3,3) line exhibits maser emission. All of the inversion 
states have similar excitation temperatures at a given density. ULIRGs such as Arp220, 
however, could show substantially different excitation temperatures due to their high gas 
densities, high kinetic temperatur e s, an d strong velocity gradients. Nevertheless, assuming 



that the IWalmsley fc Ungerechts! (119831 ) calculations provide at least an indication of the 
possible range, we can derive ammonia column densities using Eq.HJ The results are listed 
in Table[3j The individual inversion states exhibit column densities of order ~ 10 16 cm -2 for 
T ex = 50 K and ~ 10 15 cm" 2 for T ex = 10 K, which add up to (7.8 ± 0.4) x 10 16 cm" 2 and 
(1.6 ± 0.1) x 10 16 cm -2 for the two excitation temperature values, respectively. The relation 
scales linearly with T ex , so even a high value of 100 K or 150 K would only modestly increase 
the column densities. The (0,0) ortho-ammonia level is degenerate and does not exhibit an 
inversion line. However, it is the lowest energy state of ammonia and thus has a significant 
population. To estimate the amount of ammonia in the (0,0) state, we extrapolate the 
^23456 — (124 ± 19) K rotational temperature in the Boltzmann diagram (Fig.HJ assuming 
an ortho-to-para-NH 3 ratio of unity) and derive (0,0) columns of (6.0 ±0.9) x 10 15 cm -2 and 
(1.2±0.2) x 10 15 cm -2 for excitation temperatures of 50 K and 10 K, respectively. Adding this 
to the (1,1) through (6,6) columns, we derive total ammonia columns of (8.4±0.5) x 10 16 cm -2 
and (1.7 ± 0.1) x 10 16 cm~ 2 for values of T ex = 50 K and 10 K, respectively (see TableEJ). 
This value neglects all higher level populations, a valid assumption given the negligible 
contributions that these will make to the entire column if the rotational temperature stays 
the same for J > 6. For example, the (7,7) inversion line is 538 K above the ground state; 
by extrapolating T rot via Eq.EJin the Boltzmann diagram (Fig.HJ), we see that the NH 3 (7,7) 
line would fall at a N/T ex that is ~ 14 times below the NH 3 (6,6) value. 

As shown in Fig. [3] and in Tabled some peak and integrated apparent optical depths 
differ between the ATCA and GBT observations. In the following, we try to derive how 
the properties of the gas have changed if the different spectra are due to a time variation 
of the gas in front of the line of sight to the continuum (the GBT data were observed 
about a year after the ATCA data). To start, we use Eq.HJ to determine N/T ex , which is 
directly proportional to the integrated apparent optical depths. According to Table[2j this 
translates to ATCA over GBT column density ratios of ~ 1.3 and 2.6 for the NH 3 (2,2) 
and (3,3) transitions, respectively, at a given T ex . The total ammonia columns depend on 
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the temperature assumed, but, as shown in Fig.Hl the GBT data points are not reliable 
enough to determine a temperature as the NH 3 (1,1) line has very different properties. If 
the (2,2) and (3,3) lines are representative of the gas phase measured by the GBT, however, 
we can estimate a rotational temperature from the slope and we derive a cold T 2 3 ~ 42 K. 
Extrapolating this to higher transitions via Eq.[5]and adding the NH3(1,1) will then provide 
~ 28 x 10 15 cm -2 for T ex = 50 and ~ 6 x 10 15 cm~ 2 for T ex = 10, or ~ 2.5 times less than 
what we estimate from the ATCA data (see Table[3]). This is the variation of the ammonia 
column that potentially may have modified the absorption profile between the ATCA and 
GBT measurements over the time of a year. 

Eventually, we c an us e the total column densities to derive ammonia volume densities. 
Downes &: Solomon! (119981 ) show that the bulk of the molecular gas in Arp 220 is found in 
a ring or disk (the light gray shaded area in Fig.[H]). In their model, the size of this ring 
or disk is about ~ 700 pc when measured from the western edge closest to the observer 
toward the eastern nucleus, the farthest continuum source against which we see the NH3 
absorption. This geometry leads to NH 3 volume densities of fy 



0.5 



and f v 



0.5 



x (3.9 ±0.2) x 10" 5 cm- 3 

x (0.8 ± 0.1) x 10" 5 cm~ 3 for the T ex = 50 K and 10 K limits, respectively (/ v is 
the volume filling factor of the dense gas). If the line of sight through the molecular gas is 
substantially smaller, i.e. if the absorption occurs mainly against the western core of Arp 220 
(see Sect. l4.4p . the volume densities increase by a factor of ~ 3. 

We can now estimate the volume densities of molecular hydrogen that is traced by the 
ammonia a bsorp tion li nes. To do so, we ha v e to a ssume a fractional abundance of ammonia. 
Ott et al.l (120051 ) and iMauersberger et al.l (120031 ) show that for a sample of nearby, non- 
ULIRG starburst galaxies, [iV(NH3)/iV(H2)] hovers around 10 -8 . If this value is applicable 



to the dense gas in Arp 220, the H 2 densities are of order fy ' x (1 — 4) x 10 3 cm~ 3 . For 
a filling factor of unity, this is similar to t he environment of the SgrB2 region near the 
Galactic Center ( Hiittemeister et al.l 1995bh an d close to the critical den s ity of ammonia 
which is ~ 2 x 10 3 cm" 3 (IRohlfs fc Wilson! 120041 ) . IScoville. Yun. fc Bryantl (119971 ) . however, 
estimate a density of ~ 2 x 10 4 cm -3 for the molecular ring in Arp 220. If the ammonia 
absorption lines trace the same gas, the volume filling factor would be of order fy ~ 0.01. 



In Sect. 14.21 we speculate that the small absorption depth of the (1,1) line relative to the 
other NH3 lines might be caused by an additional, cold layer of ammonia in Arp 220. The 
component would have to be located inside the beam but should be displaced from the lines 
of sight toward the continuum sources and thus be seen in emission. If this line has a width 
of ~ 200 km s -1 , such as to compensate some of the absorption line, the assumed strength 
can be converted to a column density in the (1,1) line of ~ 1 x 10 14 cm~ 2 . With T12 = 20 K 



(Sect. l4~2|) and using the Equation (A15) in lUngerechts. Winnewisser. &: Walmsleyl ( I1986I ). 
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this can be extrapolated to a total column of Admission ~ 2 x 10 14 cm -2 . The hypothetical 
emission component thus has a column density more than 2 orders of magnitude smaller 
than the ammonia layer that is seen in absorption against the background continuum. 



4.4. Gas Kinematics 



The ammonia linewidths, ~ 230 — 430km s 1 (NH3(1,1): 121km s 1 ), are at the up- 
per end and even wider than those fou nd in nearby 'normal' starburst galaxies, such as 
NGC253 (~ 250 km s" 1 ; lott et al.l fioosh. Maffei2 (~ 150 km s" 1 ), or IC342 (~ 50kms -1 ; 



Mauersberger et all 120031 ) when averaged over the nuclear regions. Toward Ar p220, molec- 



ular emission lines such as CO and HCO + and HNC (e.g., see the summary in iGreve et al. 

320 — 430 km s" 1 , as are highly excited CO and CS tran- 
ce;.) which 



2008 




Greve et al. 


2009: 



show linewidths of ~ 300 km s 1 for each of the nuclear co mponents. Imaging of the CO 



(J = 2 — > 1) distributio n with sub-arcsecond resolution (jScoville. Yun. fc Bryantl Il997 



Downes fc Solomon! Il998l ) provides information on the dynamical properties of the molec- 



ular gas. Toward the east of the active region that encompasses the two galactic nuclei, 
the velocity range is 5500 — 5800 km s _1 and toward the west it is 5200 — 5500 km s~ x 



(jScoville. Yun. fc Bryantl 119971 ). One possibility is that the molecular gas we observe in am- 
monia, ~ 5450 km s _1 with a linewidth of ~ 230 — 430 km s _1 , has components toward 
both nuclei. This, however, is unlikely since the absorption spectra are well defined by a 
one-component Gaussian that cannot be explained by two disks, separated in velocity space 
well above our velocity resolution. Alternatively, the width may reflect the differenc e in 



systemic velocities of the cores, which is ~ 300 km s _1 (e.g. iDownes fc Solomon! 1 19981 ). It 
is a coincidence that the linewidths of the gas around the individual cores and the velocity 
difference between the cores is of the same order. This makes it difficult to determine the 
origin of the gas based on kinematic information alone. The NH3 spectra show similarities 
to previously measured OH lines. The absorption lines of rotationally excite d OH exhibit 
linewidths of ~ 250 - 370 km s" 1 iHenkel et al.lll986l . Il987l : ISalter et al.lboosh which is like 
the linewidths obtained from NH3 and, as NH3, OH absorption has no clear velocity imprint 
of either nucleus. Both ammonia and OH may therefore trace the same molecular material 
either from a single nucleus, such as the linewidth suggests, or, alternatively, in the common 
disk in which the two nuclei are embedded, as suggested by the systemic velocities. 

The gas in the dense cores has volume densities of ~ 10 6 cm~ 3 (see Sect. l4.3p . but 
we estimate a 2-3 orders of magnitude lower density for the NH 3 traced gas (Table[3]). A 
combination of the systemic line velocities, the derived densities, and the temperature argues 
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that the ammonia- (and OH-) traced gas is more likely to be stored in the inner ring or disk 
that surrounds both nuclei (see Fig. [6]). Unfortunately, as for ammonia, only single-dish data 
of rotationally excited OH have been published to date and its distribution with respect 
to the nuclei in Arp 220 has not yet been directly imaged. Sub-arcsecond mapping of the 
components with the EVLA or ALMA is needed to at least determine the absorption profiles 
against each nucleus. Such observations will be crucial to obtain a better model of the layered 
gas distribution in Arp 220. 



5. Summary 

In this article, we have presented ATCA and GBT observations of ammonia toward the 
ultraluminous far-infrared galaxy Arp 220. We conclude the following: 

1. We detect all observed ammonia lines ( J,K)— (1,1), (2,2), (3,3), (4,4), (5,5), and (6,6), 
toward the bright continuum in absorption. 

2. The detected ammonia linewidths are in the range ~ 120 — 430 km s _1 , in agreement 
wi th those of other mo lecules. We cannot confirm the extremely wide lines reported 



bv lTakano et al.l fl2005[ ). 



3. Based on the ATCA ammonia data (but excluding the NH 3 (1,1) line), we derive the 
rotational temperature of the molecular gas to be (124 ± 19) K. This translates to a 
kinetic temperature of T^ in = (186 ± 55) K. The kinetic temperature of NH 3 is similar 
to that found in the nucleated starburst galaxy NGC 253. 

4. We see no evidence for an ortho-to-para-ammonia ratio other than unity. This implies 
ammonia formation temperatures exceeding ~ 30 K. 

5. Combined column densities of the measured lines plus the extrapolated (0,0) column 
add up to (1.7 ± 0.1) x 10 16 cm~ 2 and (8.4 ± 0.5) x 10 16 cm~ 2 for assumed excitation 
temperatures of 10 K and 50 K, respectively. We estimate the volume density of the 
molecular gas that is traced by ammonia to be ~ /y ' 5 x (1 — 4) x 10 3 cm -3 (fy: volume 
filling factor). Together with the temperature estimate and the fact that the spectra 
are well described by single-component Gaussians, we conclude that the bulk of this 
moderately dense gas likely surrounds both starburst nuclei. 

6. Widespread gas, based on low excitation CO, exhibits densities that are about two 
orders of magnitude less than the ammonia-traced gas layer. This component defines 
the outer disk of Arp 220. In contrast, the gas close to the two starburst nuclei is 
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characterized by densities two orders of magnitude larger than the ammonia-traced 
gas. The gas that we see via ammonia is thus likely a layer that is located between 
the nuclei and the outer envelope. The geometry is most likely a ring or disk that 
surrounds both nuclei. The temperature profile appears to have a local maximum at 
the position of the gas that is traced by ammonia. 

7. The NH3(1,1) line is weaker than expected for thermal gas based on all other mea- 
sured NH3 transitions. A possible explanation for this unusual behavior could be an 
additional, cold, < 20 K component of molecular gas, located inside the beam but dis- 
placed from the lines of sight toward the continuum sources in Arp220. This gas, if 
present, would have to be at an estimated column density of A^ emission ~ 2 x 10 14 cm -2 
to sufficiently mask the (1,1) absorption. The cool temperature is required to keep the 
masking of the other lines to a minimum. 

8. We also detect an absorption feature adjacent in frequency to the NH3(3,3) line in the 
ATCA data but not in the GBT data. This absorption feature may be interpreted as 
an OH 2 n 3 / 2 J = 9/2F = 4-4 transition, or the combination of the F=4-4 and F=5-5 
OH doublet. Both options yield line parameters similar to that of ammonia, including 
a rotational temperature that is about 2 times larger for OH. If this detection holds, it 
would be the first extragalactic evidence for this line which has so far only been seen 
toward Galactic compact H II regions such as W3(OH). 

9. The ATCA and GBT data show some differences in their line shapes. This is most likely 
an instrumental effect of the two very different telescopes and observing procedures. We 
cannot exclude, however, that there is intrinsic time variability of the gas distribution 
in Arp 220 along the narrow lines of sight toward the absorbing background structure. 
If this is the case, then the GBT data suggests gas with a ~ 3 times cooler rotational 
temperature than that of the ATCA measurement, and the column density would be 
~ 2.5 times lower at the time when the GBT observations were taken (a year after the 
ATCA observations). Given the typical cloud and continuum size scales of a few to 
tens of parsec and the velocities of a few hundred km s -1 , however, drastic variations 
of the absorption spectra on annual timescales seem to be unlikely and an instrumental 
effect is more plausible. 

We thank Fabian Walter for his comments on the manuscript. C. H. thanks ATNF 
for support during his time spent at ATNF and ATCA. This research has made use of the 
NASA/IPAC Extragalactic Database (NED), which is maintained by the Jet Propulsion Lab- 
oratory, Caltech, under contract with the National Aeronautics and Space Administration 
(NASA) and NASA's Astrophysical Data System Abstract Service (ADS). 
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Fig. 1. — ATCA broadband spectrum in the 23.2 to 23.6 GHz range. Units are in apparent 
optical depth r. The bars at the bottom indicate the bandwidths and tunings of the different 
spectral settings. The rms noise of the observations is r rms O.Of 3. 
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Fig. 2. — Continuum subtracted ammonia absorption ATCA spectra of Arp220 in units of 
apparent optical depth r. Gaussian fits to the lines are overlaid. All panels are on the same 
scale. 
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Fig. 3. — GBT broadband spectrum (thick lines) overlaid over the ATCA spectrum. 
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Fig. 4. — Boltzmann plot of the ammonia measurements. The ordinate exhibits the weighted 
logarithm of N/T ex . We plot the rotational temperature fitted to all but the (1,1) datum 
with the errors in slope. ATCA data points are marked by bold circles, whereas the GBT 
data is plotted with thin crosses and errorbars. 
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Fig. 5. — LVG results on the conversion between rotational and kinetic ammonia tempera- 
tures (solid curves: para-ammonia, dotted curves: ortho-ammonia). The subscripts denote 
the transitions used for the T rot determination. The dashed straight line indicates where 
both temperatures would be identical; T rot is always underestimating Tkm- For Tyi and T 36 
we show logarithmic fits to the curves in dashed lines. These fits are described in eqs.El and 
[3 respectively. 
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Fig. 6. — A cartoon image of Arp220. Shown are the different gas phases and some charac- 
teristic molecular or atomic tracers. The outer envelope has a lower limit to the temperature 
and the NH 3 component likely resides between the outer envelope and the inner dense regions 
that surround the eastern and western cores. 
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Table 1. Observational Parameters. All observations were directed toward 
RA(J2000) = 15 /l 34 m 58!10, DEC(J2000)=23°30 / ll'. / 0. For the assumed 75Mpc distance of 

Arp220, 1" corresponds to 0.36 kpc. 



Telescope 


Date 


Center Frequencies 


BW a 




Beam size 


NH 3 Transition 






[GHz] 


[MHz] 


[MHz] 


["] 






ATCA 


2005 Jul 28-31, Aug 01 


23.14, 23.20, 23.26, 23.32 
23.38, 23.44, 23.50, 23.56 


128 


4.4 


31.1 x 21.1 


(1,1), (2,2), (3,3) 


ATCA 


2005 Jul 30 


23.665, 23.737 


128 


4.4 


31.2 x 21.4 


(4,4) 


ATCA 


2005 Jul 31 


24.052, 24.124 


128 


4.4 


30.9 x 19.8 


(5,5) 


ATCA 


2005 Mar 17 and 20 


24.603 


128 


4.4 


10.7 x 7.7 


(6,6) 


GBT 


2006 Jan 21 and Feb 02 


23.445 


800 


1.56 


33 


(1,1), (2,2), (3,3) 



a BW: bandwidth of each frequency setting 
b Au: spectral resolution 
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Table 2. Ammonia (1,1) through (6,6) line parameters. Velocities are given with respect 
to the redshift of Arp220 of z — 0.018126, or a velocity of 5434 km s _1 . The parameters of 
the OH 2 n 3 / 2 J = 9/2 F = A — 4 line candidate are also provided. 



NH 3 line 


^-peak 


peak 
opt,bary 

[km s _1 ] 


[km s _1 ] 


f rdv. 
[km s _1 ] 


ATCA 


(1,1) 


0.046 ±0.012 


63 ± 16 


121 ±38 


5.9 ±2.5 


(2,2) 


0.178 ±0.008 


-39 ±5 


232 ± 13 


43.9 ±3.2 


(3,3) 


0.145 ±0.005 


-1 ±8 


435 ± 19 


67.0 ±3.7 


(4,4) 


0.085 ±0.006 


-18 ± 10 


285 ± 25 


25.9 ±2.9 


(5,5) 


0.072 ±0.008 


45 ± 18 


327 ± 45 


25.1 ±4.4 


(6,6) 


0.53 ±0.016 


20 ±55 


268 ± 56 


15.1 ±5.6 


OH F=4-4 


0.086 ±0.006 


-13 ± 10 


292 ± 29 


26.7 ±3.2 


GBT 


(1,1) 


0.042 ±0.005 


27 ±7 


121 ± 17 


5.4 ± 1.0 


(2,2) 


0.136 ±0.004 


-35 ± 10 


230 ± 22 


33.3 ±3.3 


(3,3) 


0.094 ±0.003 


-42 ± 10 


258 ± 22 


25.8 ±2.4 
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Table 3. Densities of the molecular gas. The center column is derived using an excitation 
temperature of 50 K and the right column using T ex = 10 K. Whereas Ej=i 6 N(J, J) 
denotes the sum of the measured lines, the ammonia column densities listed as 
Ej=o 6 J) include the (0,0) inversion state that was extrapolated using T rot = 124 K. 
The column densities transform into volume densities nNH 3; j=o...6 for an assumed line of 

sight of 700 pc. Finally, for a fractional abundance of [A^(NH 3 )/A^(H 2 )] = 10~ 8 , the 
resulting molecular hydrogen densities % 2 are shown in the last row. fy is the volume 

filling factor. 



T ex = 50 K T ex = 10 K 



N [10 15 cm- 2 ] 



NH 3 (1,1) 


3.7 ± 1.6 


0.7 ±0.3 


NH 3 (2,2) 


21.0 ± 1.5 


4.2 ±0.3 


NH 3 (3,3) 


28.3 ± 1.6 


5.6 ±0.3 


NH 3 (4,4) 


10.1 ± 1.1 


2.0 ±0.2 


NH 3 (5,5) 


9.3 ±1.6 


1.9 ±0.3 


NH 3 (6,6) 


5.3 ±2.0 


1.1 ±0.4 


NR 3 (0,0) extr 


6.0 ±0.9 


1.2 ±0.2 


Ej=i... 6 N(J, J) 


78 ±4 


15.5 ±0.9 


Ej=0...6 N ( J i J ) 


84 ±5 


16.7 ± 1.0 



/ v - 5 x n N H 3 ;J=o...6 [10 5 cm 3 ] 
3.9 ±0.2 0.8 ±0.1 



/ v a5 xn H2 [10 3 cm- 3 ] 
3.9 ±0.2 0.8 ±0.1 



